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Abstract 

We present a calculation of e + e~ production in 7 A reactions at MAMI and 
TJNAF energies within a semi-classical BUU transport model. Dilepton in- 
variant mass spectra for 7C, ^yCa and 7PD are calculated at 0.8, 1.5 and 2.2 
GeV. We focus on observable effects of medium modifications of the p and u 
mesons. The in- medium widths of these mesons are taken into account in a 
dynamical, consistent way. We discuss the transport theoretical treatment of 
broad resonances. 
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I. INTRODUCTION 



In-medium properties of hadrons are of fundamental interest with respect to an under- 
standing of QCD in the non-perturbative regime (cf. []]]). During the past decade especially 
the properties of vector mesons have found widespread attention as they may be related to 
chiral symmetry 

A comparison of experimental data on dilepton production in nucleus-nucleus collisions 
at SPS energies with transport theoretical calculations [J7L|S| seems to indicate a lowering 
of the p-meson mass in the nuclear medium. However, since in a heavy ion collision the final 
dilepton yield is obtained by an integration over dileptons emitted at different densities and 
temperatures a discrimination between different scenarios of in-medium modifications for 
the vector mesons is difficult PJlQ"l. 

Moreover, there is a more fundamental concern: ultrarelativistic heavy ion reactions pro- 
ceed, at least in their initial stages, far from equilibrium, whereas all theoretical predictions 
of in-medium properties are based on equilibrium assumptions. Therefore it is necessary to 
probe in-medium properties of vector mesons under 'cleaner' conditions. For that purpose 
photon or pion induced reactions are promising tools. In such reactions the nuclear medium 
is very close to equilibrium at normal nuclear density and temperature zero. The predicted 
in-medium effects for the vector mesons are so large that they should have observable con- 
sequences already for densities p = po- Here one should note that even the dileptons seen in 
ultrarelativistic heavy ion reactions stem to a large part from densities p < 2po- 

Our calculations are based on a semi-classical BUU transport model that has recently 
been very successfully applied to the description of heavy ion collisions at SIS energies 
| lT|l and photoproduction of pions and etas in nuclei fT2|. Meanwhile we have extended 
the model to the description of dilepton as well as strangeness production and to the high 
energy regime by including the Lund string fragmentation model FRITIOF ||13|| . This allows 
us to calculate inclusive particle production in heavy ion collisions from 200 AMeV to 200 
AGeV, in photon and in pion induced reactions with the very same physical input. To 
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our opinion, the simultaneous description of as many experimental observables as possible 
is necessary because of the large number of parameters, like unknown cross sections, and 
the strong assumptions that enter semi-classical transport models. In heavy-ion collisions 
particle production depends not only on a large number of elementary reaction channels 
but also on the global space-time evolution of the system. On the other hand photon, 
pion or proton induced reactions allow to fix certain ingredients of the transport model 



as observables depend in general only on a few reaction channels. In Ref. |14| we have, 
for example, determined the in-medium 7/-nucleon cross section from photoproduction of 
?7-mesons. Just recently we were able to improve considerably our treatment of the A- 
resonance by comparing our calculations to experimental data on photoproduction of pions 

Within our model we have already given predictions for dilepton production in pion 
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nucleus reactions [16j that will be measured by the HADES collaboration at GSI fL7| . In 
the present paper we want to look at e + e~ production in photonuclear reactions in the 
energy range from 0.8 to 2.2 GeV that will be accessible at TJNAF [T8|] and at its lowest 
energy also at MAMI. 

Our paper is organized as follows: In Section II we describe our model. Here we focus 
on the treatment of broad resonances and our description of elementary photon nucleon 
reactions. In Section III we present our results for dilepton production in photon nucleus 
reactions and discuss the possibility to subtract the Bethe-Heitler contribution. In Section 
IV we include in-medium modifications for the vector mesons into our calculations and 
present their effect on the dilepton yield. We close with a summary in Section V. 

II. THE BUU MODEL 

The transport model used here has been developed starting from the model that has 



been described in full detail in Refs. [11, IS]. Here we restrict ourselves to the description of 



the essential new features of our method. 



3 



A. Resonance properties 



Instead of the baryonic resonances in Ref. [ I L| that were taken with their parameters 
from the PDG |2(J we now use the resonances from the analysis of Manley and Saleski |[21|| . 
This has the advantage of supplying us with a consistent set of resonances. In particu- 
lar, now the resonance parameters are consistent with the parameterizations of the reso- 
nance widths. We take into account all resonances that are rated with at least 2 stars 
in Ref. 0: P 33 (1232), Pn(1440), P>i 3 (1520), 5 U (1535), P 33 (1600), S 3 i(1620), 5 n (1650), 
D 1B (1675), F 1B (1680), Pi 3 (1879), S 31 (1900), P 35 (1905), P 31 (1910), P 35 (1930), P 37 (1950), 
Pi 7 (1990), Gi 7 (2190), P 35 (23 50). The resonances couple to the following channels: Nit, 
Nr} : Nuj, AK, A(1232)tt, Np, Na, A^(1440)tt, A(1232)p. The cross section for the produc- 
tion of a resonance R in a collision of a meson m with a baryon B is given by: 

2J R + 1 47T ar^rff 
mB ^ R (2J m + 1)(2J B + 1) P ( a - MlY + s rgf 2 ' 1 ) 

where Jr, J m , Jb denote the spins of the resonance, the baryon and the meson, respectively. 
k is the cms momentum of the incoming particles, s is the squared invariant energy, and 
Mr is the pole mass of the resonance. The total decay width T°" t * is given as sum over the 
partial decay widths of the resonance. For a specific channel mB it is: 

-pout _ -pO PmB\S) / s 
mB mB i -\r \ j l z J 

with being the partial decay width at the pole of the resonance and p m B{s) is given as: 

PmB{s) = J dp m dp B A m {p m )A B {pB) q ^ S '^^ B ^ B^ mB (qR) , (3) 

where A m and Ab denote the spectral functions of the outgoing particles, q is their cms 
momentum, l m B is their relative orbital angular momentum, and B\ mB is a Blatt-Weisskopf 
barrier penetration factor p2"| , plf (interaction radius R = lfm). For the spectral function 
Ai of an unstable particle i we use: 

A(u) = - VlIjeM (4) 
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where Mi denotes the pole mass and r to4 (/z) is the total width. Here we neglect any spin 
degrees of freedom as well as a momentum dependence of the real part of the self-energy. 
For a stable particle (with respect to the strong interaction) we simply have: 

A i (jj) = 6(ji-MJ . (5) 

The incoming width T™ B in Eq. ([!]) is given by: 

T n kBf (kR) 

-ptn riiR pO hnB V > (c\ 

1 mB — ^mB L mB /- / fl/f \ ) W 

y sp mB \M R ) 

where is the appropriate Clebsch-Gordan coefficient for the coupling of the isospins of 
the baryon and the meson to the isospin In of the resonance. 

We note that we use - in contrast to Ref. |21| - relativistic propagators in Eqs. ([[]), and 
momentum dependent widths in the spectral functions of the outgoing particles. However, 
this has only a very small effect on the resonance production cross sections and does not 
require a readjustment of the resonance parameters. 

The mesonic resonances are treated analogously to the baryonic ones, i.e. their two-body 
decay widths are calculated according to Eq. (fj). The used parameters are given in Table |. 
The three pion decay width of the tu-meson is assumed to be constant since it is very small. 

B. The collision term 

1. Baryon-baryon collisions 



For invariant energies < 2.6 GeV we describe baryon-baryon collisions as in Ref. [11 
with the same matrix elements. Our modified treatment of the resonance properties pre- 
serves the very good agreement with the experimental data on one- and two-pion production 



in nucleon-nucleon collisions shown in Ref. ||11|| . For higher energies we use the string frag 



mentation model FRITIOF |H| with the same parameters as in Ref. [p3 |. This approach is 



similar to the ones in the hadronic transport models described in Refs. [p4| , |25| . 
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2. Meson-baryon collisions 



For meson-baryon collisions we use the string fragmentation model for s/s > 2.2 GeV. For 
lower energies the most important contributions come from intermediate nucleon resonances 
which are described according to Eq. (|l]). 

In case of 7r~p-scattering the incoherent sum of all resonance contributions gives a very 
good agreement with experimental data for pion momenta p n < 1.1 GeV, corresponding to 
invariant energies \/s < 1.73 GeV. In Fig. [I] we show the total n~p cross section. For higher 
energies (1.73 GeV < yfs < 2.2 GeV) we additionally include a background ttN — > 7r7riV cross 
section in order to reproduce the total cross section for which we use a parameterization 
from the PDG (shown as solid line in Fig. |]). In Fig. ||] we show that the resonance 
contributions give a very good description of experimental data on elastic scattering, charge 
exchange, two-pion and eta production cross sections. 

For 7r + p-scattering one gets, also for lower energies, a good agreement with experimental 
data only if one takes into account all resonances from Ref . [^TJ , including the 1 star ones, as 
can be seen from Figs. |] and [|. Since we do not explicitly propagate the 1 star resonances 
we put their contributions into background cross sections for ttN — > ttN and itN — > tckN. 
For higher energies we, again, include a two-pion production background term that is fitted 



to the total cross section from Ref. 

For 7r°p-scattering we have from isospin symmetry: 

1 

for the total cross section. The cross sections for pion-neutron scattering also follow from 
isospin symmetry. 

In Fig. |^ we compare the resonance contributions to 7r~p — > np° to the experimental data 
from Ref. ||27|| . While for invariant energies \fs above 1.8 GeV the experimental data are 
reasonably well described by the resonances there is a strong disagreement at lower energies, 
in particular in the region of the _Di3(1520)-resonance. The p-mesons produced at these 
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energies have invariant masses essentially below the pole mass m° p . In Ref. |27 the p-meson 
cross section has been obtained by a fit to invariant mass spectra of the outgoing pions in 
7i~p — > n7T + 7T~ . For low \fs the shapes of the p- meson and the 'background' contributions 
become very similar and a determination of the p-meson cross section gets very difficult. 

In Ref. |21]] the couplings of the baryonic resonances to the Np channel have been de- 
termined by using amplitudes for ttN — > Np that were obtained by a partial wave analysis 
of all available ttN — > Ntcit data in Ref. ||28|| . We also note that the large coupling of 



the .Di3(1520)-resonance to the Np channel found in Ref. [21] is in line with other similar 
analyses ||29|| . Therefore we consider the experimental data in Ref. [p7| to be wrong for low 

In addition to the resonance contributions we include the following processes: 



tcN «- 


-> luN 


nN - 


^ ujttN 


uoN - 


-> ttttN 


uN - 


-> uN 


ttN «- 


-> 4>N 


vriV - 


-> (pirN 


<j)N - 




<j)N - 





where we adopt the cross sections from Ref. |30[. For cross sections involving an to- meson 
we, of course, subtract our resonance contributions from these cross sections. 



C. Treatment of broad resonances 



In our model broad resonances, like the baryon resonances or the p-meson, are not just 
produced and propagated with their pole mass but according to their spectral function. The 
transport equation for a system of N particle species reads: 



where Fi(f,p, p,, t) denotes the one-particle spectral phase space density of particle species i 
with r and p being the spatial and momentum coordinates of the particle. p is the invariant 
mass of the particle and Hi(f,p, p, F\, . . . , Fj?) stands for the single particle mean field 



Hamilton function which, in our numerical realization [ IT] , is given as: 



Hi= y/QjL + Stf + p* , (8) 

where Si(f,p,p,Fi,...,F N ) is a scalar potential. We note that we neglect a vec- 
tor potential and Si is an effective scalar potential that is - for the nucleons 
obtained from a non-relativistic potential (for details see Ref. [ITp. The terms 



Gi(f,p, p,, Fi, . . . , F N ) and Li(r,p, p, F 1; . . . , F N ) stand for a gain and a loss term, respec- 
tively, and Ai(f,p, p, Fi, . . . , Fn) is the spectral function of particle i. The distribution 
function fi is defined by: 

fi{r,p,fj.,t) = — - ; 9 
Ai{r,p,n,t) 

for stable particles it reduces to the usual phase space density. 

In order to be more specific about Gi and Li let us consider, as an example, a system 
of nucleons, rho-mesons, pions and a single baryonic resonance species R that are coupled 
only via R «-> Np and p — > tttt. Then the gain term G p is: 

^^i/^^^^'^'^^^^-^'^'^ ' (10) 

where the factor (1 — f n ) accounts for the Pauli principle of the outgoing nucleon and 
Pn = Pr — P P due to momentum conservation. For simplicity, we have neglected a possible 
finite width of the nucleons as well as a Bose enhancement factor (1 + f p ) for the p-meson. In 
the differential decay width dT R _ Np the spectral function of the p-meson A p enters simply 
as a multiplicative factor (Eqs. (0),(0)). Therefore, here G p does not depend on A p since we 
neglect any process where more than one p-meson is produced, like e.g. R — > ppN . 
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The loss term L p reads: 




r 




3 fnyj Vni t)^np^pn—^R i 



(11) 



with v np being the relative velocity of the nucleon and the p-meson and <r pn ^R their cross 
section for the production of resonance R (Eq. ([[])). Tp^^ denotes the two-pion decay width 
of the p-meson in the calculational frame. 

The total in-medium width T to t, P appearing in the spectral function from Eq. @) is 
directly related to the loss rate L p : 



where 7 is a Lorentz factor which appears since T tot p is the decay rate in the rest frame of 
the p-meson. 

The loss and gain term for the resonance R can be written down in an analogous way. 
One immediately sees that the transport equations of the p-meson and the resonance R 
are coupled in a higly non-linear way. Especially the in-medium widths, that are functions 
of space time and 4-momentum, of both particles depend on each other through integral 



The above described equations can easily be extended to a transport model with more 
particle species and a realistic collision term. We also note that it is straightforward to 
formulate the theory consistently for real and imaginary part of the self-energies of the par- 
ticles. However, in our model we treat real and imaginary part of the self-energies completely 
independently. This violates analyticity, but in the present stage it would already require a 
considerable effort to treat the imaginary parts of all particles in a realistic transport model 
in a completely self-consistent way. 

In Ref. [|3T] it has recently been stressed that, because of unitarity, it is important to 
respect Eq. ( |T2"D in transport calculations. This means that in the population of a resonance 
the same width has to be used in the spectral function that enters the dynamical calculation 
via the collision term. In Refs. [0,^] we have already taken into account the in-medium 



(12) 



equations (Eqs. © , © , © , (|) , (0) ) . 
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widths of the baryonic resonances A(1232), iV(1520), iV(1535), and iV(1680) in a consistent 
way for population and propagation. Since during a photoproduction reaction the nucleus 
remains, in the time interval relevant for meson production, close to its ground state the 
calculational effort is in this case manageable. As we describe the nuclear ground state in 
a local Thomas-Fermi approximation we can use nuclear matter values for the in-medium 
widths that depend only on the invariant mass p, the absolute value of the 3-momentum 
and the density p: 

r(f,t,p,n) -»• r(p(?,t), \p\,n) . 

In the present paper we only take into account the in-medium widths of the p- and 
cu-mesons. In particular, we neglect any medium-modifications of the iVp-widths of the 
baryonic resonances that would, in a self-consistent calculation of the self-energies, directly 
follow from a modification of the spectral function of the p-meson. As was shown in Ref. |33| 



such effects might give large modifications of the baryon widths and also influence the 
spectral function of the p-meson. However, an inclusion of such effects would enhance the 
numerical effort dramatically, especially if one takes also the modifications of the real parts 
of the self-energies into account. Moreover, we do not expect such effects to be relevant for 
photon energies above 1 GeV because here the nucleon resonances that lie below + m° p 
and might thus get strongly modified, like the _D 13 (1520), play only a minor role. Observable 
effects of different medium modifications of the D ri will be reported elsewhere. 

The transport equations Eq. (0) do not yet give the correct asymptotic spectral phase 
space densities for particles that are stable in vacuum. This can be seen by noting that 
a collision broadened particle does not automatically lose its collisional width when being 
propagated out of the nuclear environment. The same problem appears for resonances whose 
imaginary part of the self-energy is non-zero in-medium in kinematical regimes where it is 
zero in vacuum. 

The reason for this deficiency is directly related to the semi-classical approximations 
on which the transport equation Eq. fl7|) is based, in particular the neglect of coherence 
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effects. However, since a realistic quantum transport theory is numerically not yet realizable 
we will nevertheless work with Eq. (|7]) as it is certainly a step beyond the usual on-shell 
approximation. Moreover, if the time evolution of the system is such that the rate of particle 
production and absorption is much larger than the change of the width with time the problem 
with surviving off-shell contributions will be negligible. Under the assumption that the gain 
term in the collision term is in magnitude comparable to the loss term we can formulate this 
in the following way: 

rdt»^ . (13) 

In case of a particle moving in a static nuclear medium with density profile p(r) we can 
rewrite this condition: 

(M) 

P A 

where e p is a unit vector along the momentum direction of the particle and A denotes its mean 
free path. In Section |1V A| we will discuss the validity of this condition for our calculations 



and present possible remedies to the problem. 

D. Parameterization of the elementary jN cross sections 

For invariant energies \/s < 2.1 GeV, corresponding to E 1 < 1.88 GeV on a free nucleon 
at rest, we describe one-pion, two-pion and eta production as in Ref. Hl2| . For the two- 



pion production cross sections on the neutron we meanwhile use the experimental data from 



Refs. |3^J35[1 instead of the recipe described in Ref. fl2|| . For larger energies we use, like 
for the hadronic interactions, the string fragmentation model FRITIOF where we initialize 
a zero mass p°-meson for the photon following a VMD picture. For the total cross section 



we use a parameterization from Ref. [Q. The Lund model is then used to determine the 
probabilities for the different final states. In Fig. |^ we show that this gives a very good 
description of charged particle multiplicities in photon-proton collisions; the agreement seen 
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there is better than could be expected from a model that had been developed for applications 
at high energies. However, we do not expect the Lund model to give correct predictions for 
all specific channels, especially with respect to isospin. The role of the Lund model for our 
calculations is to supply us with an overall description of the elementary reaction dynamics 
in order to allow us to take into account multi-step processes where, for example, a primary 
produced pion produces a vector meson on a second nucleon. 

The vector meson (p,u,<f>) production in 7iV — > VN collisions is fitted to experimental 
data f36|,|37}] and treated independent of the Lund model also for high energies. The cross 



section is given as: 



— f dfi\M v \ 2 PfA v (fi) , (15) 



(7 7 jV->jVV 

V 

where y/s is the total energy of the jN system, Pi,pf are the momenta of the initial and 
final particles in the center-of-mass system, and Ay is the spectral function of vector meson 
V (Eq. U). The matrix elements Aiy are parameterized in the following way: 

|.M p | 2 = 0.16 mb-GeV 2 (16) 

\M^\ 2 = 0.004 mb-GeV 2 . 

In Fig. [7| we show the resulting cross sections (dash-dotted curves) for jp — > pp° (upper part), 
7P — > pu (middle part) and 'jp —>■ p<p (lower part) in comparison with the experimental data. 
For the angular distribution of the produced vector mesons we use 

dcr , , 

— ocexp(^) , (17) 

where t denotes the square of the 4-momentum transfer of the photon to the vector meson. 
In Ref. ]36[ the parameter B was, dependent on photon energy, fitted to p°-production. 
Here, we adopt these values and also use them for uj- and ^-production. 

In our calculations there is an additional contribution to 7iV — > Np° coming from the 
decays of the iV(1520) and iV(1680) resonances which is shown in Fig. [7] by the dashed line. 
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These decays predominantly contribute to low mass p-mesons below the experimentally seen 
p production threshold. 

Besides the exclusive process — > VN we also have, for the photon energies considered 
here, to take into account additional channels for photoproduction of vector mesons. For 
energies above 2.1 GeV we calculate these cross sections using the Lund model. Below 2.1 
GeV we absorb everything into the channel 7iV — > VttN for which we use the following cross 
section: 

16(2tt) 7 



' 'yN — *VttN 



■ J dpd$ 3 \M V7r \ 2 A v (p) , (18) 



where d^ 3 denotes the 3-body phase space element as, for example, given by Eq. (35.11) 



in Ref. [20]. The matrix elements are adjusted to give a continous transition to the string 



fragmentation model at \fs = 2.1 GeV. We use: 

\M p o n \ 2 = IM^l 2 = 0.5 mb . 

The calculated inclusive rho and omega production cross sections (after the subtraction of 
the exclusive part) are shown as dotted lines in Fig. [7|. The total inclusive vector meson 
production cross sections are indicated as solid lines. 

For photon energies above 1 GeV we take into account nuclear shadowing of the incoming 
photon by adopting the model of Ref. ||38|| . Since this has, for the photon energies considered 
here, only a small impact on the final results our specific transport theoretical realization 



will be described elsewhere 39 . 



E. Dilepton production 



In our analysis we calculate dilepton production by taking into account the contributions 
from the Dalitz-decays A — > Ne + e~ , r] —>■ / ye + e~ , to — > ir e + e~ , tt° — > 7e + e~ and the direct 
dilepton decays of the vector mesons p,u>,(fi. 

The Dalitz decays of the 7r° and the i] are parameterized according to Ref. [[TI|]. For the 



Dalitz decay of the to we use the parameterization from Ref. ||40|| . The A Dalitz-decay is 
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described in line with Ref. ||2| where we, however, use g = 5.44 for the coupling constant 



in order to reproduce the photonic decay width r (0) = 0.72 MeV. 

The dilepton decay of the vector mesons is calculated assuming strict vector meson 
dominance |E| as in Ref. 0: 
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77? 

r v ^ e+e -(M)=C VJ f 3 , (19) 
where C p = 8.814 x 10~ 6 , C u = 0.767 x 10" 6 and C = 1.344 x 10" 6 , respectively @. Within 



an extended vector meson dominance picture |44| one has a dilepton decay amplitude that 



consists of 2 terms, one describing the coupling of the virtual photon to the hadron with a 
strength proportional to M 2 and another with strength proportional to M°. However, since 
we neglect a direct coupling of the virtual photon the use of strict vector meson dominance is 
more appropriate. We also note that this dilepton decay width together with our parameters 
for the p-meson gives a very good description of the experimental data for e + e _ — > ?t + ?t~. 

In our transport model the dilepton yield is obtained from the phase space distributions 
of the respective sources by a time integration. For the vector mesons the massdifferential 
dilepton production is given as: 

%^=r^ r -^F y (f, ( , P »^ , (20) 
a/J Jo (2vt j d 7 

where 7 is a Lorentz factor which appears since Ty^ e + e - is the width in the rest frame of 
the vector meson. The Dalitz decay contributions contain an additional mass integration. 
For the A-resonance we have for example: 

dN A ^ Ne+e - f°° ^ 3 d 3 p 1 dT A ^ Ne+e - 



> dt d r ——-r dfi 2 F A (r } t } p,iJ,2) . (21) 

a ix Jo (^n) 6 7 dfx 

III. DILEPTON PRODUCTION IN jA REACTIONS 

In Figs. |5| [TTJ| we present the calculated dilepton spectra da/dM for 7C, 7Ca, 7Pb reac- 
tions at photon energies E 1 = 0.8, 1.5, 2.2 GeV. A mass resolution of 10 MeV is included 
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through a convolution of our calculated spectrum with a Gaussian. Here neither collisional 
broadening nor a mass shift of the vector mesons were taken into account. 

The thin lines indicate the individual contributions from the different production chan- 
nels; i.e. starting from low M: Dalitz decay 7r° — > / ye + e~ (short-doted line), rj — ► r ye + e~ 
(dotted line), A — > Ne + e~ (dot-dashed line), uj — > 7i°e + e~ (dot-dot-dashed line); for M ~ 
0.8 GeV: p° — > e + e~ (dashed line), u — > e + e~ (dot-dot-dashed line), — > e + e~ (dashed 
line). The full solid lines represent the sum of all sources. The dominant processes in the 
low mass region up to M ~ 500 MeV are the T], u and A Dalitz decays. Above M ~ 0.6 GeV 
the spectrum is dominated by the vector meson decays with a low background from other 
hadronic sources. 

In our calculations we only take into account p-mesons with masses above which is 
the threshold of the strong decay since in the calculation of the spectral function (Eq. (|4])) 
we neglect contributions from electroweak decays. Therefore we get a discontinuity of our 
spectra in Figs. |8HT0| at the two-pion mass that is, however, because of the other sources 



and the mass resolution hardly visible. 

A. Bethe-Heitler contribution 

Besides the 'hadronic' contributions as discussed above we also have to take into account 
dilepton production via the so-called Bethe-Heitler (BH) process for which the Feynman 



diagrams are depicted in Fig. |TT| that contribute to lowest order in the electromagnetic 
coupling constant a. 

On a single nucleon, with the electromagnetic form factors known from electron scatter- 
ing, the BH process is completely determined by QED dynamics and can easily be calculated. 
For a detailed description of the involved matrix elements we refer to Ref. |PJ from which 
we also adopted the parameterizations of the electromagnetic form factors Wi(Q 2 ,u) and 
^(Q 2 , v) of the nucleon. 

For our calculations we take only the incoherent sum over BH contributions on single 
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nucleons into account and neglect contributions where the intermediate photon couples to the 
charge of the whole nucleus. While the latter will, because of the Z 2 dependence, dominate 
all integrated cross sections it can experimentally easily be suppressed by appropriate missing 
mass cuts. 

In Fig. [T2| we compare the BH contributions for 7Pb at 1.5 (upper part) and 2.2 GeV 
(lower part) with the 'hadronic' contributions that we have already shown in Fig. 0. One 
sees that the sum of the elementary cross sections on the proton and the neutron (dashed 
lines) is much larger than the 'hadronic' contributions (solid lines) for dilepton masses below 
0.6 GeV. In the region of the p and uj meson the BH contribution is about a factor of 4 smaller 
but still non-negligible. With the inclusion of Fermi motion and Pauli blocking (dotted lines) 
the BH contribution is reduced significantly for low invariant masses but hardly affected for 
masses larger than 700 MeV. 

In order to suppress the BH contribution we implemented the cuts (k ■ p),(k ■ p+) > 
0.01 GeV 2 , where k is the 4-momentum of the incoming photon, and p, p + are the 4-momenta 
of electron and positron, respectively. These cuts, which reflect the pole-like behaviour of 
the intermediate electron propagator, suppress the BH contribution by a factor of 10 - 
dot-dashed lines in Fig. O - and practically do not have any influence on the 'hadronic' 
contributions. 

In our calculations we do not take into account interference terms between the BH and 
the 'hadronic' contributions. For an inclusive cross section that contains to each e + e~ pair 
the pair with exchanged momenta the interference term vanishes and the BH contribution 
can easily be subtracted. Therefore we will in the following discuss only the hadronic 
contributions. 



IV. IN-MEDIUM EFFECTS IN DILEPTON PRODUCTION 
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A. Collisional broadening 



The in-medium widths of the p and uj mesons are calculated as sketched in Section [11 C| . 
In the rest frame of the meson the total in-medium width is given as: 

rf ot Gu, \p\, P ) = v v vac {p) + r^fa \p\, P ) , (22) 

where the collisional width T^ oll reads: 

rl M (ji,\p\,p)='yp<v V N(T$ t N> , (23) 

and T^ ac is the vacuum decay width. In Eq. (p3|) the brackets stand for an average over the 
Fermi sea of the nucleons, vvn is the relative velocity between vector meson and nucleon, 
and ay l N is their total cross section, p is the nuclear density and 7 the Lorentz factor for 
the boost to the rest frame of the vector meson. 

The upper part of Fig. ^ shows the collisional width of the p meson as a function of 
momentum and mass at nuclear matter density p = po- The structure at low p comes from 
the resonance contributions, especially from the Di3(1520). Note that the width becomes 
very large (up to about 600 MeV), corresponding to a complete melting of the p-meson. 
The lower part of Fig. [L3| shows the momentum and density dependence of the uj collisional 
width calculated at the pole mass p = . At nuclear matter density po an d a momentum 
of p = 1 GeV we obtain a collisional width of about 80 MeV which is about a factor of 10 
larger than the vacuum decay width. 

In Fig. [14| (upper part) we show the contribution of the p-meson to the e + e~-yield for 
7Pb at 1.5 GeV. The solid line indicates the bare mass case (as in Fig. |H]), i.e. without 
collisional broadening. The curve labelled 'coll. broadening' (dash-dotted line) results if we 
take into account the collision broadening effect in the production of the p-mesons. Here we 
calculate the p-meson production cross sections in photon-nucleon collisions (Eqs. (|I5|) ,( |I8|) ) 
with the in-medium spectral function. Since the in-medium spectral function depends on 
the momentum of the p-meson with respect to the nuclear medium an additional angular 
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integration has to be performed. For the case of exclusive production *yN — > Np we use 
the angular distribution from Eq. (|ITD . For — > Npn we assume an isotropic three-body 
phase space distribution. 

As discussed in Section [11 C| we do not modifiy the iVp-width of the baryonic resonances 
but the masses of the p-mesons stemming from these decays are distributed according to 
the (phase space weighted) in-medium spectral function. We, again, neglect p-mesons with 
masses below the two-pion threshold. 

From Fig. [14] one sees that the inclusion of collisional broadening leads to a depletion 
of the p-meson peak by about 30% and a shift of strength to lower dilepton masses. One 
also observes a very large peak at M — 2m^ which is in fact a pole but here finite due to 
our numerical solution. The reason for this divergence is directly related to our discussion 
in Section [II C| of the asymptotic solutions of the semi-classical transport equation when 
including in-medium spectral functions. At the two-pion mass the vacuum spectral function 
of the p-meson is zero while the in-medium spectral function has some finite value since the 
collision width from Eq. ( P3|) does not vanish. When travelling to the vacuum the respective 
component of the spectral phase space density becomes infinitely long lived and leads to 
the divergence. If we included the electroweak decay width of the p-meson into the collision 
term of the transport equation the problem would not be solved. The pole would only be 
replaced by a numerically indistinguishable large peak. 

In Fig. [14] (lower part) we show the contribution of the tu-meson when including collisional 
broadening (dash-dotted line) in comparison to the calculation with the vacuum spectral 
function (solid line, as in Fig. [10]). One observes a strong broadening of the u;-peak which is 
partly covered up by the inclusion of a mass resolution of 10 MeV. However, such a strong 
broadening is not realistic because most of the tu-mesons that contribute to the dilepton 
spectrum decay in the vacuum. This is also reflected by the violation of the condition 



Eq. ( |i~4"|) for the validity of our off-shell transport equation in this case. 
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1. Prescriptions to obtain correct asymptotic behaviour 

In the following we want to discuss a prescription that allow us to obtain a divergence free 
p-meson and a reasonable u;-meson contribution. For that purpose we introduce a potential 
that shifts a particle to its vacuum spectral function when it propagates to the vacuum. 
Such a potential can not be defined on the level of the transport equation Eq. (|7|). However, 
we can introduce such a potential on the level of our numerical realization. We recall that 
we solve the transport equation by a so-called test particle method, i.e. we make an ansatz 
for the spectral phase space density F: 

F(F,*,^^)ocX;^-^(*)W-ft(*))^-^) , (24) 

i 

where fi(t), Pi(t), and /ij denote the spatial coordinate, the momentum and the mass of the 
test particle i, respectively. Now, we can define for each test particle a density dependent 
scalar potential Sj in the following way: 

s;(Pi(0) = (/C ed -/C c )^ , (25) 

Pi 

where pJ™ ed is the 'in-medium' mass of the test particle chosen according to the mass dif- 
ferential production cross section. p^ ac is the 'vacuum' test particle mass which is chosen 
according to the production cross section with a vacuum spectral function. p^ r is the baryon 
density at the creation point, whereas Pi(t) = p(r*j(t)) is the baryon density during the prop- 
agation. Setting the test particle mass \x\ to the 'vacuum' mass p\ ac the effective in-medium 
mass p* is given as: 

tiMt)) = + 8i (p(t)) . (26) 

Eq. (p6|) gives the correct asymptotic behaviour: the effective mass at the creation point 
is kept unchanged (/x*(pf ) = p™^)', during the propagation the test particle mass changes 
linearly with density and outside the nucleus it becomes equal to the bare mass p*(pi = 
0) = pl ac . The potential s« enters into the test particle propagation as a usual potential and 
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therefore guarantees that this prescription does not violate energy conservation. [] 

The potential Eq. ( f25|) can assume rather large values in the case of broad resonances 
even if the in-medium corrections are small. Its effect is, nevertheless, in this case negligible. 
As a check we have performed a calculation of photoproduction of p-mesons in which we 
distributed the masses yJ^ ed according to the vacuum spectral function. This gave practically 
the same result as the calculation without potential since the lifetime of the p-mesons is so 
short that only very few propagate through a relevant density gradient for which the potential 
becomes important. 



In Fig. [L4| we show the results of our calculations with the described prescription for 
the p and the u mesons (curves labelled 'propagation to bare mass', dashed lines). For the 
p-meson the divergence at the two-pion threshold is removed. For invariant masses above 
500 MeV we get practically the same result as without the prescription because of the reason 
described in the preceding paragraph. The broadening of the cu-peak is reduced because now 
only the cu-mesons that decay inside the nucleus contribute to the broadening. 

Another possibility to overcome the problems with the transport theoretical treatment 
of broad resonances is to avoid their explicit propagation. Therefore we have also performed 
simulations in which we calculated the cross sections for elementary dilepton production via 
vector mesons as instantaneous one-step processes. The results are shown in Fig. [14| by the 
curves labelled 'instantaneous decay' (dot-dot-dashed lines). The p-meson contribution is 
almost identical to the one calculated with the prescription described above because of its 
very short life time. For the a;- meson we get a reduction by about a factor of 3 since in 
the instantaneous decay scheme we neglect the possibility that an tu-meson can escape from 
the nucleus and then contribute to dilepton production with a much larger branching ratio 
than inside the nucleus. As such dynamical effects are important to take into account we 



1 After this paper was submitted two other papers dealing with the propagation of broad resonances 
appeared on the preprint server [46,47]]. 
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consider the description via an instantaneous process to be inadequate. 

The easiest way to cope with the divergence of the p-meson contribution is to use a 
minimum two-pion decay width. In Fig. [TJ] we show the result of a calculation where we 
set this minimum width to 10 MeV (dotted line). One sees that still a large peak remains. 
Moreover, this approach is in any case questionable. 

For the reasons described above we consider the prescription 'propagation to bare mass' 
as the only possibilty to implement collision broadening effects in our calculations. However, 
we want to stress that this prescription is not fully satisfactory since it is only formulated 
on the level of our specific numerical solution of the transport equation and since it neglects 
any memory effects. 



B. 'Dropping' vector meson mass 

In order to explore the observable consequences of vector meson mass shifts at finite 
nuclear density the in- medium vector meson masses are modeled according to Brown/Rho 
scaling [|J or Hatsuda and Lee f|] by introducing a scalar potential Sy(r): 

S v (r) = -am v ^ , (27) 
Po 

where p(f) is the nuclear density, my the pole mass of the vector meson, p = 0.168 fm~ 3 , 
and a ~ 0.18 for the p and u. The effective mass p* is then given as: 

fi* = n + S v . (28) 

For the effective pole mass Pq we thus get: 

In photon-nucleus reactions vector mesons are produced with large momenta relative 
to the nuclear medium. Within a resonance-hole model for the p-meson self energy in the 
nuclear medium it has been shown in Ref. E8[ that the real part of the in-medium p-meson 



self energy increases with momentum and crosses zero for a momentum of about 1 GeV. In 
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order to explore the implications of such a behaviour we also use a momentum dependent 
scalar potential S™ om for the vector mesons: 



In our calculations we take the vector meson potentials into account for the calculation of the 
phase space factors in — > NV (Eq. (|15D ) and — > NVtt (Eq. (|18D ). We neglect these 
modifications for the vector meson production in the string fragmentation model FRITIOF 
and also do not modify the iVp-widths of the baryonic resonances. 



In Fig. [T5| we show the contribution coming from the a;-meson to e + e~-production in 7Pb 
at a photon energy of 1.5 GeV. A dropping mass scenario according to Eq. fl2~7|) (dot-dashed 
line) gives a two peak structure corresponding to tu-mesons decaying inside and outside 
the nucleus, respectively. An additional inclusion of collisional broadening, as described in 



inside the nucleus. The height of the peak around the vacuum pole mass of the cu-meson 
is hardly affected by the dropping mass scenario. On the one hand the lowering of the 
mass reduces the vacuum peak because tu-mesons decaying inside the nucleus contribute to 
lower masses. On the other hand the total production of u;-mesons is enhanced since the 
phase space factors entering the elementary photoproduction cross sections (Eqs. (|15"D, ([T8|)) 
are increased for lower masses. In our calculation both effects nearly cancel each other for 
masses around the vacuum pole. 

In Fig. [16] (upper part) we show the total e + e~ yield for the same reaction. A dropping 
mass scenario for the vector mesons (dashed line) leads to a second peak structure at in- 
variant masse of about 650 MeV. The peak around 780 MeV remains practically unchanged 
since it is dominated by cu-mesons decaying outside the nucleus. 

With the inclusion of collisional broadening the in-medium peak gets completely washed 
out (dotted line). The dilepton yield at intermediate masses is about a factor of 2 larger 




(29) 



C. Dileptons from reactions: In-medium effects 



Section [IV A| gives a substantial broadening of the dilepton yield from cu-mesons that decay 
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than in the bare mass case. At the two-pion threshold there is a visible discontinuity which 
results from our neglect of p-mesons with effective masses below the two-pion mass. 

Using the momentum dependent scalar potential from Eq. ( p9|) we obtain the curves 
labelled 'momentum dependent potential' (dot-dot-dashed lines). The result is very close to 
the bare mass case as the vector mesons are mainly produced with momenta around 1 GeV 
for which the potential is zero. 

In the lower part of Fig. [16| we show that the effect of a dropping vector meson mass at 
a photon energy of 2.2 GeV is qualitatively the same as for 1.5 GeV. The calculation with a 
momentum dependent potential gives again a result that is practically the same as for the 
bare mass case. 

In photonuclear reactions vector mesons are in general produced with larger momenta 
relative to the nuclear medium than in heavy-ion collisions. Since the in-medium spectral 
functions of the vector mesons are momentum dependent one might thus observe rather 
different in-medium effects in both reactions. These, together with additional information 
from pion- nucleus and proton- nucleus reactions |T7|, might help to discriminate between 
different scenarios of medium modifications. Therefore a calculation of all reactions within 
one model is necessary for a conclusive interpretation of the experimental data. Our BUU 
transport model provides such a tool. 



V. SUMMARY 

We have studied e + e~ production in 7C, 7Ca, and 7Pb reactions at photon energies of 
0.8, 1.5, and 2.2 GeV within a semi-classical transport model. Various contributions were 
taken into account for dilepton production: Dalitz decays of A(1232), 7r°, rj, and u as well 
as direct dilepton decays of the vectors mesons p, u, and <p. We have focused on observable 
effects of in-medium modifications of the vector mesons p and u. 

It was shown that the Bethe-Heitler process which dominates all integrated cross sections 
for dilepton production can be sufficiently suppressed by appropriate cuts on the lepton 
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momenta. For dilepton invariant masses above 600 MeV the spectrum is dominated by the 
the decays of the vector mesons (p, lu, (p). A mass shift of these mesons as proposed in 
Refs. |3],|4[] leads to a substantial enhancement of the dilepton yield at invariant masses of 
about 650 MeV by about a factor of 3 and should clearly be visible in experiments that will be 



carried out at T JNAF |T3] . However, a calculation for which we used a linearly momentum 



dependent potential for the vector mesons |48j gave practically no effect compared to the 
bare mass case. 

We have stressed the necessity of a simultaneous description of vector meson produc- 
tion in different nuclear reactions as one probes in-medium properties at different momenta 
relative to the nuclear medium. 

Exemplarily for the case of the p and the u mesons we have discussed the conceptual 
problems in the treatment of broad resonances in semi-classical transport models. We have 
presented a prescription that allows to obtain reasonable results when taking into account 
in-medium spectral functions. 
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TABLES 

TABLE I. Properties of mesonic resonances. Mq and Tq denote the pole mass and the width 
at the pole mass, respectively. 



Meson 


M [MeV] 


L [MeV] 


decay channels 


P 


770 


151 


TTTT 




782 


8.4 


TTTT (2%), 7T°7 (9%), TT+TT-TT (89%) 




1020 


4.4 


P TT (13%), KK (84%), tt+tt-tt (3%) 


a 


800 


800 


TTTT 



28 



FIGURES 




FIG. 1. The total it p cross section compared to the experimental data from 
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FIG. 2. 7r p cross sections for elastic scattering, charge exchange, two-pion and eta production. 



The experimental data are taken from [49]. 
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FIG. 3. The total ir + p cross section compared to the experimental data from fl49| . 
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FIG. 4. 7r + p cross sections for elastic scattering and two-pion production. The experimental 
data are taken from 
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FIG. 6. Charged particle multiplicity cross sections in 7p reactions: jp — > 1 charged particle 
(dashed line (calculation as described in the text), squares (experimental data from Ref. |4J|)), 3 
charged particles (dotted line, circles), 5 charged particles (dot-dashed line, up triangles), 7 charged 
particles (dot-dot-dashed line, down triangles). Also shown is the total cross section (solid line, 
rhombs). 
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FIG. 7. The p (upper part), u (middle part) and <j) (lower part) meson production cross sec- 
tion for 7p reactions. The experimental data corresponding to the exclusive reactions "fp — > Vp 
(V 



p,uj,(j)) are taken from Ref. [36] (full circles) and from Ref. |37j (open squares). The 
dash-dotted lines are our parametrization of the exclusive data; the dashed line (upper part) is the 
resonance contribution. The dotted lines indicate the calculated inclusive vector meson produc- 
tion cross section(see text); the solid lines correspond to the total vector meson production cross 
section. 
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FIG. 8. The dilepton invariant mass spectra da/dM for 7C at the energy of E 1 = 0.8 GeV 
(upper part), 1.5 GeV (middle part) and 2.2 GeV (lower part) calculated with bare meson masses 
including a mass resolution of 10 MeV. The thin lines indicate the individual contributions from 
the different production channels; i.e. starting from low M: Dalitz decay 7r° — > je + e~ (short-doted 
line), r\ — > 7e + e~ (dotted line), A — > Ne + e~ (dot-dashed line), oj — ► 7r°e + e _ (dot-dot-dashed line); 
for M « 0.8 GeV: w — > e + e~ (dot-dot-dashed line), p° — > e + e~ (dashed line). The full solid line 
represents the sum of all sources. 
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FIG. 9. The dilepton invariant mass spectra da/dM for ^yCa at the energy of E~ = 0.8 GeV 
(upper part), 1.5 GeV (middle part) and 2.2 GeV (lower part) calculated with bare meson masses 
including a mass resolution of 10 MeV. The assignment is the same as in Fig. || 



35 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 

M [GeV] 



FIG. 10. The dilepton invariant mass spectra da/dM for 7PD at the energy of £L = 0.8 GeV 
(upper part), 1.5 GeV (middle part) and 2.2 GeV (lower part) calculated with bare meson masses 
including a mass resolution of 10 MeV. The assignment is the same as in Fig. |8|. 
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FIG. 12. The dilepton invariant mass spectra da/dM for 7PD at the energy of E~ = 1.5 GeV 
(upper part) and 2.2 GeV (lower part). The solid lines indicate the 'hadronic' spectra as in 
Fig. [l(^. The dashed lines are the Bethe-Heitler contribution calculated as a sum of incoherent 
contributions from 7+nucleons. The dotted lines show the BH yield with taking into account 
nucleon Fermi motion and Pauli blocking. The dot-dashed lines are the BH terms with the cuts 
(k-p),(k-p+) >0.01 GeV. 
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FIG. 13. The upper part shows the collisional width of the p meson as a function of momen- 
tum and mass at normal nuclear density p = po- The lower part is the momentum and density 
dependence of the oj collisional width calculated at p = m°. 
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FIG. 14. The dilepton yield from p (upper part) and oj (lower part) meson decays for 7PI) at 
1.5 GeV calculated within different prescriptions: the solid lines are the results with bare mass, the 
dot-dashed lines indicate the calculation with collisional broadening, the dot-dot-dashed lines cor- 
respond to the instantaneous decay, the dotted line (upper part) is the result with T^ in = 10 MeV, 
the short dashed lines indicate the yield calculated with the 'propagation to bare mass' prescription. 
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FIG. 15. The dilepton yield from u mesons for 7PD at 1.5 GeV. The solid line indicates the 
bare mass case, the dot-dashed line is the result with in-medium masses, the dashed line shows the 
effect of collisional broadening together with the dropping mass. 
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FIG. 16. The dilepton invariant mass spectra for 7PD at 1.5 GeV (upper part) and 2.2 GeV 
(lower part). The solid lines indicate the bare mass case, the dashed lines are the result with 
the dropping mass scenario, the dotted line (upper part) shows the effect of collisional broadening 
together with the dropping mass, the dot-dot-dashed lines indicate the result with the momentum 
dependent potential from Eq. (pS| ). 
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